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Abstract 

Shell model calculations are performed for magnetic dipole excita¬ 
tions in ®Be and ^*^Be in which all valence configurations plus 2hL0 ex¬ 
citations are allowed (large space). We study both the orbital and spin 
excitations. The results are compared with the ‘valence space only’ 
calculations (small space). The cumulative energy weighted sums are 
calculated and compared for the J = O"'' T=0 to J = 1'^ T=1 exci¬ 
tations in ®Be and for J = O'*' T=1 to both J = 1+ T=1 and J=l'^ 

T=2 excitations in ^'^Be. We find for the J = O"'' r=l to J = 

T=1 isovector spin transitions in ^*^Be that the summed strength in 
the large space is less than in the small space. We find that the high 
energy energy-weighted isovector orbital strength is smaller than the 
low energy strength for transitions in which the isospin is changed, 
but for J = 0"*“ T=1 to J = T*" r=l in ^*^Be the high energy strength 
is larger. We find that the low lying orbital strength in ^*^Be is anoma¬ 
lously small, when an attempt is made to correlate it with the B{E2) 
strength to the lowest 2"*“ states. On the other hand a sum rule of 
Zheng and Zamick which concerns the total B{E2) strength is rea- 





sonably satisfied in both ®Be and ^^Be. The Wigner supermultiplet 
scheme is a useful guide in analyzing shell model results. In and 
with a Q ■ Q interaction the T = 1 and T = 2 scissors modes are 
degenerate, with the latter carrying | of the T = 1 strength. 
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1 The Experimental Situation 


From our perspective, much experimental information is lacking in the 
nuclei ®Be and ^°Be. For example, no J = states have been identihed in 
^°Be. The B{E2) from the 2i state of ®Be to the J = O’*" ground state is not 
known -this is understandable because of the large decay width to two alpha 
particles. 

The following states and their properties are of interest to us: 

(a) ®Be 

The J = 2f state has an excitation energy of 3.04 MeV. The J = Af state 
is at 11.4 MeV. This is consistent with an J( J-l- 1) spectrum of a rotational 
band, but it should be recalled that any spin-independent interaction gives 
an J(J -|- 1) spectrum in the p shell. The J = 1^ T = 1 state, which we 
discussed extensively in a previous publication is at 17.64 MeV and the 
J = l^" T = 0 state is at 18.15 MeV. 

The B{M1) from the 17.64 MeV state to the ground state has a strength 
of 0.15 W.u. or B{M1)\= 0.27/rAr^. The B{M1) of this state to the 2f state 
is 0.12 W.u. or B{M1)\= 0.21pN^ [|]. Of course B(M1) |=3i?(Ml) | 

(b) i°Be 

The J = 2 1 state is at 3.368 MeV and the J = 2^ state at 5.960 MeV . 
We recall that with a spin independent interaction the 2f and 2^ would be 
degenerate. The experimental spectrum looks more vibrational. However, 
the values of B{E2)\iTom the J = 0^ ground state to the 2f state is very 
strong: B{E2)W 52 e^fm^. Raman et. ah deduce from this a deformation 
parameter (3 = 1.13 |^. As mentioned above, there are no J = I"*" states 
mentioned in the compilation of Ajzenberg-Selove [Q. Also the J = 4+ state 
has not been found. 
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2 The Interactions 


We have chosen two types of interactions to do the calculations. First we 
use a short range ‘simplified realistic’ (x, y) interaction previously used by 
Zheng and Zamick |Q, and then we use a long-range quadrupole-quadrupole 
interaction. By choosing these two extremes, we make sure that the results 
we obtain are not too dependent on the specifics of the model. 

In more detail, the(a;, y) Hamiltonian is: 

n = T.T.+T.ntj) 

i<j 

where 

V = V, + xVso + yVt 

with c =central, s.o. =spin-orbit, and t =tensor. 

For (x, |/) = (1,1) the matrix elements of this interaction are close to those 
of realistic G matrices such as Bonn A. We can study the effects the spin-orbit 
and tensor interactions by varying x and y. 

Note that we do not add any single-particle energies to the above Hamil¬ 
tonian. Rather, we let the single-particle energies be implicitly generated by 
H. Hence, if we set x=0 i.e. turn off the two-body spin-orbit interaction, 
we will also be turning off the one-body spin-orbit splitting coming from this 
interaction. 

As a counterpoint, we repeat all the calculations with the Q ■ Q Hamilto¬ 
nian 

= E ^ -xT.Q-Q 

i i<j 

Note that we have added the term which is not present for the (x, y) 

interaction. The reason for this is that Q ■ Q cannot generate any single¬ 
particle potential energy splitting whereas the (x, y) interaction can. 


4 



Whereas the {x, y) interaction like all realistic interactions is of short 
range, the Q ■ Q interaction is long range. Yet, as we shall see some of the 
results (but not all) are rather similar for the two interactions. Since the 
best milieu for the existence of scissors mode excitations (orbital magnetic 
dipole excitations) are strongly deformed systems, one would expect the Q ■ Q 
Hamiltonian to yield strong scissors modes. But is this also true for the 
realistic interaction ? We will address this question. Another motivation for 
introducing the Q ■ Q Hamiltonian is that it is easy to establish a connection 
via energy weighted sum rule techniques between isovector orbital H(Ml)’s 
and isoscalar and isovector B{E2)^s. 

We shall be performing the calculations, not only in the Oku space (small 
space) but also in a space which allows 2hu admixtures (large space). For 
the Q ■ Q Hamiltonian in the small space the energy matrix is proportional to 
y. Hence the energy eigenvalues depend linearly on y, but the eigenfunctions 
(and H(Ml)’s and B{E2)^s) are independent of the interaction strength. In 
a large space calculation there is one more parameter; the energy splitting 
induced by ^ i.e. 2hu. Thus the wave function and the corre¬ 

sponding i?(Ml)’s and B{E2)^s will also depend on y. 

We have chosen values of y appropriate for the large space calculation. 
We also use these same values in the small space. One can argue that in the 
small space one should use a renormalized value y' which is close to twice 
y. However, as mentioned above, the wave function and hence B{M1) and 
B{E2) will not change, only the energies. By choosing the same y in the two 
spaces it is easier to see what the differences in the two calculations are. The 
values of y are 0.5762 for ^Be and 0.3615 for ^^Be. 


5 



3 The summed magnetic dipole strength 


In Table I we give the summed magnetic dipole strength {J2i B{M1 : 
0+ T = 1 ^ 1+ T = 1) and Ei B{M1 : 0+, T = 1 ^ 1+, T = 2) ) 
broken up into isoscalar and isovector and spin and orbit and where we use 
the {x,y) interaction with x = 1, y = 1. We hrst discuss the behaviour as 
a function of the size of the model space. Later we will make a comparison 
of the behaviour in ®Be and ^°Be. There are striking differences for the two 
nuclei. 

Our small space calculation consists of all conhgurations of the form 
(0s)"^(0j9)^ for ®Be and (0s)"^(0j9)® for ^'^Be. The large space consists of those 
conhgurations plus 2hu excitations. Thus one can either excite two particles 
to the next major shell or excite one particle through two major shells. We 
also give results for the summed strength in the low-large space -this is the 
low energy part of the large space covering an energy range more or less equal 
to that of the small space. It is easy to identify the low energy sector because 
there is a fairly wide plateau in the summed strength which separates the 
low energy rise from the high energy rise. 

Usually the large space summed strength is somewhat larger than the 
small space strength e.g. for the isovector orbital strength in ®Be the values 
shown in Table / are 0.6701 and 0.7283 respectively. But there is 
one glaring exception. For the case of = 0^ T = 1 —= 1+ T = 1 
transitions in ^'^Be, the summed isovector spin strength in the large space is 
2.08 10“^ ylf but in the small space it is bigger 2.34 10“^ fj,%. For the orbital 
strength it is the other way around but for the physical case {giT, = 1, g^, = 0, 
gsTT = 5.586, gsu = —3.826) the spin prevails and the summed strength in the 
large space 1.952 is less than in the small space 2.09 

Thus it is not always true that the net result of higher shell admixtures 
is to rob strength from the low energy sector and move it to higher energies. 
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In some cases the total strength gets depleted. 

We next compare the low energy sum in the large space with the small 
space sum. In all cases the latter is larger than the low energy sum, thus 
indicating that there is a quenching of the low energy part due to higher 
shell admixtures. The hindrance factor [{low large)/small] is 0.88 for the 
isovector orbital in ®Be, 0.73 for the isovector spin in ®Be, 0.77 for the total 
Ml in ^°Be etc. 

Note that the total Ml’s for ^°Be are somewhat larger than for ®Be. How¬ 
ever there is a dramatic drop in the orbital strength in ^°Be relative to that 
in ®Be. The large space summed orbital strength for ®Be is 0.73 whereas 
for ^°Be (to J = 1+ T = 1 and T = 2) the value is (0.196 -|- 0.183)=0.38 
/i^. In the low energy sector the ®Be value is 0.59 gi\ whereas for ^°Be it is 
0.10 -|- 0.13 = 0.23 less than half the value for ®Be. 

From the systematics of orbital transitions in heavy nuclei one concludes 
that the proper milieu for isovector orbital transitions is strongly deformed 
nuclei. Can one conclude that ^°Be is not strongly deformed? The answer, 
by examining the tables of Raman et. ah is no! There is a strong E2 
connecting the 0^ and 2f states in ^°Be. From this the authors conclude 
that the deformation parameter (3 is about 1.13 -quite enormous. Of course 
®Be might have an even stronger E2 transition -there is no data on this in 
the Raman paper [Q, probably because of the rapid decay of the 2f state 
into two alpha particles. 

4 The cumulative energy weighted strength 
for orbital transitions in ^Be and ^^Be 

In this section we present results and hgures for the cumulative energy 
weighted sum of magnetic dipole strength. 
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We are motivated in so doing by various energy-weighted sum rules that 
have been developed e.g. by Zheng and Zamick p, Heyde and de Coster 
[^], Moya de Guerra and Zamick [^, Nojarov |P, Hamamoto and Nazarewicz 
[^] and Fayache and Zamick |^. We will focus in particular on the orbital 
strength for which the operator is (L^r — Ly)/2. In a previous publication 
we presented results for the {x,y) interaction with x=l, y=l for ®Be |]^. In 
this work the quadrupole-quadrupole interaction results are compared with 
the {x, y) interaction results, and furthermore we extend the calculation to 
^°Be. In the latter nucleus one does not have N = Z and this leads to big 
differences. 

Whereas in ®Be there is only one isospin channel for isovector transitions 
J = O^" T = 0 —J = !■'■ T = 1, in ^°Be there are two: J = 0^^ T = 1 — 

J = 1+ T = 1 and J = 05^T = 1—i>J = l+T = 2. The low lying J = 1+ 
T = 1 states in ^°Be are expected to have much smaller excitation energies 
than the J = !■*■ T = 1 states in ®Be. This makes it easier to look for such 
states experimentally. 

In Table II we present the results for the summed energy weighted strengths 
for the (x, y) interaction. As a crude orientation it should be noted that sim¬ 
ple models e.g. the Nilsson model used by de Guerra and Zamick [0 and a 
model by Nojarov [§] would have the ‘large’ result be twice the ‘low large’ 
result. On the other hand Hamamoto and Nazarewicz have argued that 
the ‘large’ result should be much more than twice the ‘low large’ result. The 
actual ratios for the (x, y) and Q ■ Q interactions for this calculation (all Oho; 
conhgurations plus 2h,oj excitations) are 


RATIO 


^Be J = 0+ T = 0 ^ 
J = 0+ T = 1 ^ 
J = 0+ T = 1 ^ 

lo^e (T = 1’ + 






{x,y) 

Q-Q 

J = 

1+ 

T = 

= 1 

1.75 

1.37 

J = 

1+ 

T = 

= 2 

2.00 

1.52 

J = 

1+ 

T = 

= 1 

3.22 

3.68 

T = 

:2’) 


2.52 

2.33 


For we should actually compare the theoretical models with the com¬ 
bined result (‘T = 1’ + ‘T = 2’). 

These results indicate that the simple models are not too bad as a first 
orientation but there are fluctuations -sometimes the ratio is less than two, 
sometimes greater. We will discuss these matters in more detail in the context 
of the figures. 

In Figs 1 and 2 we show the cumulative energy weighted isovector orbital 
strength distributions in ®Be for the (x, y) interaction and for the Hamiltonian 
Hq i.e. the quadrupole-quadrupole interaction. The results for the two 
interactions are quite similar. The outstanding feature is that there are two 
rises separated by a rather wide plateau. For the {x, y) interaction the first 
rise is to a plateau at about 12 yj^MeV followed by a second rise to about 
20.8 y?pfMeV. For the Q ■ Q interaction the first plateau is at 10.25 fij^MeV 
and the second at 14 y^jMeV. A simple self-consistent Nilsson model was 
shown to give the second plateau at twice the energy of the first plateau 0 
[^]. That is to say the high energy rise was equal to the low energy rise. 
In the more detailed calculations performed here the high energy rise is less 
than the low energy rise. 

We next turn to ^^^Be. Here there are two channels: J = 0^ T = 1 — 
J = 1+ T = 1 and J = {i'^T = l^J=V^T = 2. Let us discuss the 
latter channel first. The behaviour for T = 2 in ^‘^Be is similar to that for 
T = 1 in ®Be. As shown in Figs 3 and 4 for the (x, y) and Q ■ Q interactions 
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respectively, there are two rises separated by a plateau and here the second 
rise is about twice the hrst rise for the {x, y) interaction. For the Q ■ Q 
interaction (with y = 0.3615) the low energy rise is to 1.7 y%MeV and the 
next rise is to 2.6 yj^MeV -only 1.5 to one. 

There is a big difference in the cumulative energy weighted distributions, 
shown in Figs 5 and 6, for the J = 0’''T=1—>J=1’''T = 1 channel. For 
the {x, y) interaction the hrst plateau (at about 2.5 yj^MeV) is not very hat, 
but the most outstanding feature in the curve is that the high energy rise 
is much larger than the low energy rise. The energy weighted sum reaches 
up to about 8 y%MeV. Thus the high energy rise is over three time the low 
energy rise. For the Q ■ Q interaction, the hrst plateau is better dehned -it is 
located at 1 yj^MeV and the cumulative sum extends to about 3.8 y^MeV. 


5 The Zheng-Zamick Sum Rule 


Energy weighted sum rules for magnetic dipole transitions, be they spin 
or orbital, are highly model dependent. An energy weighted sum rule for 
isovector orbital magnetic dipole transitions for the quadrupole-quadrupole 
interaction Q-Q was developed by Zheng and Zamick p. This was motivated 
by the work of the Darmstadt group showing a linear relationship 

between summed orbital B{M1) strength and the square of the deformation 
parameter i.e. 5^. 

The result was 


Y,{En-E„)B{Ml% = ) 0i isoscalar -^(-^2, 0i ^ ‘^i)isovector] 


where B{Ml)o is the value for the isovector orbital Ml operator {gi-„ = 
0.5 giu = —0.5 gs-K = 0 gsv = 0) and the operator for the E2 transitions 


{EWSR) 
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is J 2 protons epr‘^Y2 + 'Eneutrons enr‘^Y2 witli Cp = 1, 6^ = 1 for the isoscalai 
transition, and Cp = 1, = — 1 for the isovector transition. 

Let ns now describe in detail how this snm rnle works. The snm rnle 
shonld work for single-shell as well mixed-shell space. 

We hrst consider the case of ®Be. We have the following valnes in a 
large space calcnlation for the Hq interaction corresponding to orbital Ml 
excitations from the J = O’*" T = 0 gronnd state to all 1+ T = 1 states: 


1. Energy weighted isovector orbital Ml strength: 

2. The isoscalar snmmed strength B{E2; 1,1): 

3. The isovector snmmed strength B{E2; 1,-1): 

4. The right hand side = 0.1032): 


14.040 

237.46 e^fm^ 
89.611 e^fm^ 
15.25 n%MeV 


We don’t get exact agreement (14.04 n%MeV vs. 15.25 fi^MeV) bnt it 
is reasonably close. One possible reason for the disagreement is that spnrions 
states have been removed and/or that only 2huj excitations to the AN = 2 
shell are taken into acconnt p 

There have been other approaches, especially in the context of the Inter¬ 
acting Boson Model [^] which relate the energy weighted orbital magnetic 
snm to the B{E2) of the lowest 2~^ state. As a matter of cnriosity we shall 
examine in onr calcnlation what happens if we take only the lowest 2~^ state 
in the right hand side of the snm rnle {EWSR). 

For the 2/ state in ®Be we obtain (in onr large space calcnlation) B{E2] 1,1) 
196.76 and B{E2-,1, —1) = 0 (becanse the 2^ state has T = 0). The right 
hand side becomes 20.30 fi%MeV. We get a larger answer nsing the lowest 
2+ state than we do if we nse all 2+ states in the Oku and 2hoo space. The 
reason for this is that when the lowest 2+ state is exclnded, the isovector 
B{E2) is larger than the isoscalar B{E2). 
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Can we also get agreement for the sum rule in the small space (Os)^(Op)"^ 
for ®Be? Now the numbers are: 


1. Energy weighted isovector orbital Ml strength; 

2. The isoscalar summed strength B{E2; 1,1): 

3. The isovector summed strength B{E2; 1, —1): 

4. The right hand side = 0.1032); 


6.411 jjj^MeV 
72.54 

10.37 e'^fm^ 
6.414 /i^MeC 


We get perfect agreement. 

We next consider ^°Be. The relevant numbers for the large space calcu¬ 
lation are: 


1. Energy weighted Ml strength: J = O’*", T 


J = 0+, T 

Left Hand Side 


2. E5(E2;1,1) 

J = 0+, T 

3.(a) EB{E2-1,-1) 

J = 0+, T 

(h)EB{E2;l,-l) 

J = 0+, T 

(c) EB{E2-1,-1) 

Total 


4. Right hand Side = 0.0647): 


1 ^ J = 1+, T = 1 
1 ^ J = 1+, T = 2 

1 ^ J = 2+, T = 1 
1 ^ J = 2+, T = 1 
1 ^ J = 2+, T = 2 


3.811 jjj^MeV 
2.602 n%MeV 
6.413 nj^MeV 
251.4 e'^fm^ 
94.02 fm^ 
48.78 e^fm^ 
142.8 

7.029 n%MeV 


For ^“^Be we are also curious to see what happens if we use only the lowest 
2~^ state in the right hand side of the sum rule. But we have to be careful! It 
turns out that there is substantial B{E2) strength to the two lowest J = 2+ 
states. This can be understood from the fact that with a Q ■ Q interaction 
in a small space calculation ((Os)^(Op)®) the two lowest 2+ states are exactly 
degenerate. The states belong to the [/] = [42] representation. The Q ■ Q 
interaction fails to remove the degeneracy of these states. Another way of 
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stating this is that the (A/i) values for both states are (22), and the allowed 
values of the K quantum number in the Nilsson scheme are K = n, fi — 2, 
etc. Thus the 2~^ states have = 0 and K = 2. 

When we go to the large space calculation with &Q-Q interaction, limiting 
the excitations to 2huj, the degeneracy is removed but the states are still fairly 
close together. The calculated values are: 

h;2(l,l) E2{l,-l) 

2t E = 2.08 MeV 64.94 12.32 

2^ E = 2.92 MeV 93.38 10.11 


Thus, using the calculated values of B{E2) for the lowest two 2+ T = 1 
states in ^°Be, we get for the right hand side of the sum rule a value of 8.80 
H%MeV. Again, as in the case of ®Be, this is larger than the value 7.03 
H%MeV that is obtained by using all 2~^ T = 1 and all 2'^ T = 2 states. 

The corresponding numbers in small space for ^°Be are: 


1. Energy weighted Ml strength: J = O’*", T 


J = 0+, T 

Left Hand Side 


2. E5(E2;1,1) 

J = 0+, T 

3.(a) EB{E2;1,-1) 

J = 0+, T 

(b)E5(i?2;l,-l) 

J = 0+,T: 

(c) EB{E2;1,-1) 

Total 

4. Right hand Side 

= 0.0647): 


1 ^ J = 1+,T = 1 
1 ^ J = 1+, T = 2 

1 ^ J = 2+,T = 1 
1 ^ J = 2+, T = 1 
1 ^ J = 2+,T = 2 


0.7597 

H%MeV 

1.266 

^ilMeV 

2.026 

fi%MeV 

68.31 


33.80 


3.203 

e^fm^ 

37.003 

e^fm^ 

2.026 

fi%MeV 
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6 A discussion of the calculated B(E2) values 


Although the main thrust of this work is on B{M1) values, we have 
established a connection with B{E2) for the orbital case. It is therefore ap¬ 
propriate to discuss the calculated B{E2) values -comparing the behaviours 
in ®Be and ^°Be, and comparing the different interactions that have been 
used (see tables III and IV.) 

In making the comparison between ®Be and ^°Be we should lump together 
the B{E2)^s of the first two 2+ states in ^°Be because with the interactions 
used here -especially Q ■ Q- these states are nearly degenerate. (However, 
experimentally the states are well separated E 2 + = 3.368 MeV and E 2 + = 
5.958 MeV). When this is done we find that the B{E2) values in the two 
nuclei are comparable. 

For the {x,y) interaction, the calculated (large space) value of B{E2) to 
the lowest two 2+ states in ^°Be is 22.90 whereas it is 25.97 e^fm^ 

to the lowest 2’*' state in ®Be. With the Q ■ Q interaction the two values 
are respectively 46.40 and 49.16 e^frn^. One big difference between the 
two interactions is the ratio of large to small space values for corresponding 
B{E2) values. In ®Be the ratio of the large sum to the small sum is 1.98 for 
the {x, y) interaction whereas it is much larger 3.28 for the Q ■ Q interaction. 
There is much more core polarization with the Q ■ Q interaction than with 
the {x, y) interaction. 

There have been many discussions concerning the correlation of summed 
orbital Ml strength and the B{E2) from the J = 0^ ground state to the 
first 2+ state. The latter is an indication of the nuclear deformation. We 
have noted that the calculated values of B{E2) are about the same in ®Be 
and ^°Be. Thus we would expect the orbital Ml strengths in the two nuclei 
to be about the same. 

There is a certain ‘vagueness’ in what is meant by ‘strength’. It is clear 
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that the experiments thus far sample only low energy strengths up to about 
6 MeV in heavy deformed nuclei H H. Also some of the theories involve 
summed strength per se and others involve the energy weighted strength. 
Rather than enter into deep philosophical discussions about what is meant 


by strength, we will give a variety of ratios of strength 


loRe 


in Table V. We 


see that all the ratios, be they non-energy weighted or energy weighted, be 
they in small spaces or in large spaces, are substantially less than one. In 
forming the ratios, we added for the numerator (^°Be) the J = O’*" T = 1 to 
J = 1+ T = 1 and J = 0’''T = ltoJ=l+T = 2 strengths. 


7 A comparison of the J = 1+ ^ 0^ and J = 1+ 
^ 2i Magnetic Dipole Transitions 

Let us assume that the 0^ and 2f states are members of a K = 0 
rotational band and that the states have K = 1. We can then use the 
rotational formula of Bohr and Mottelson (Eq. 4-92) in their book |T^ (Ki = 

Ji J2 J 

Ml M2 M 


0, K 2 = 1) (We use the notation 
coefficient): 


for a Clebsch-Gordan 


{K2l2\\M{\)\\Ki = Oh) = ^2(2/1 + 1) 


h A I 2 
0 Ko Ko 


{K2\M{X,u = K2)\Ki = 0) 


From this we can easily deduce 

B{M1)j=i+^X=1^J=2+,K=0 


r = 


1 

2 


(Ml) j=i+^7^=1^ j=o+,/s'=o 

Note, however, that the experimental ratio for from the J = 1+ T = 1 
state at 17.64 MeV (see section 1) is = q_ 3 _ Bohr and Mottelson later 
discuss corrections to the above simple formula. 
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We can obtain a value for the above ratio by forming an intrinsic state 
and projecting out states of angular momentum J = 0 and J = 2. We 
assume an extreme prolate shape for ^Be and put the four valence nucleons 
(A^|,7V|,P|,PJ,)in the lowest Nilsson orbit with A = 0. The asymptotic 
wave function is NrYiQ. 

The J = 0 wave function is 


n'Y. 

L 


I I L 
0 0 0 


1 1 L 
0 0 0 


L L 0 
0 0 0 


[LLf 


0.74536[0 0]° + 0.66666[2 2]° 


where the notation [L^ L^Y is used. 
The J = 2 wave function is 


’ 1 1 0 ' 


’112’ 

[0 2]2 + 

’112’ 


’110’ 

0 0 0 


0 0 0 

0 0 0 


0 0 0 


[2 0]2 


’112’ 


’112’ 


- 1 

to 

to 

to 

_1 

0 0 0 


0 0 0 


000 


^•^"=2 = 0.62361[0 2]2 + 0.62361[2 0]^ + 0.47141[2 2]^ 


We don’t actually have to specify the J = 1+ state in detail to get the 
ratio. We note that only the component [2 2 ]^ of the J = 1 wave function can 
contribute to the Ml transition. The probability of this component factors 
out in the ratio. With some additional Racah algebra, we can show that 
r = I = 0.875. This should be compared with the value r = | of the simple 
rotational formula and with the experimental value r^xp = 0.8. We see that 
we get better agreement by using this projection method. 

To complete the story, using results described in the next section, we are 
able to deduce that for ®Pe P(Ml)^+^ 2 + = £;: Aw 
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8 Supermultiplet Scheme with a. Q ' Q inter¬ 
action 

8.1 A. Supermultiplet Scheme in ^Be 

The Q ■ Q interaction that we have been nsing hts in nicely with the L — S 
supermnltiplet scheme of Wigner |^. For the p shell the nnitary gronp U{3) 
is relevant since there are three states: L = 1, M =1, 0 and -1. A very usefnl 
reference for this section is the book by Hammermesh H. 

If the Hamiltonian were a Casimir operator of U{3) all states of a given 
special symmetry [/] = [/i, / 2 , /s] wonld be degenerate. For the case of 
Ip shell a state with a given particle symmetry [/i, / 2 , /s] is identical to 
a qnantnm oscillator symmetry state [^, |T^ (A,/r) = (/i — / 2,/2 — /s)- 


The states (A, p, L) are eigenstates of our Q ■ Q interaction which is a linear 
combination of the Casimir operator of SU{3) and an L ■ L interaction. The 
latter gives rise to a terminating rotational L{L + 1) spectrum for states of 
different L but with the same [/]. Amusingly, as has been pointed out by 
many, one gets identical bands in all p shell nuclei with this model provided 
the coefficient of L ■ L is hxed. 

Unlike in the s, d shell, nothing new is added by using the quantum 
numbers (A, /i) instead of [/i, / 2 , /s] for Ip shell states. This is because 
the number of different M states availbale for particles (3) coincides with the 
number of possible directions for oscillator quanta (aj,, a^, a^) and a single 
creation operator correponds to each particle. 

In more detail, the Casimir operator is C 2 = Q ■ Q — 3L ■ L. Hence, 


{~xQ ■ Q)x L — X —{C2)\ fj, + 3L{L + 1) 


= X 


-4(A^ + /i^ + A/i + 3(A + p) + 3L{L + 1) 
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where x = X§^ with 6 the harmonic oscillator length parameter (6^ = ^). 
The magnetic dipole modes in the L S T representation are: 

L = 1S' = 0T = 0 L = 0S'=1T = 1 (isovector spin mode) 

L = 1S' = 0T = 1 (scissors mode) L = 15'=1T = 0 

L = 0^ = 1T = 0 L = 1S =1T = 1 


With the Q ■ Q interaction that we have chosen, transitions from the 
L = 0 S' = 0 ground state in ^Be to all of these modes except one will 
vanish. The only surviving mode is the L = 1S' = 0T = 1 scissors mode. 

Let us give briefly the energies and some properties of the states in ^Be 
(X = 0.1865): 

(a) [f] = [4,0] (A,/i) = (4,0) Ground State Band 


L 


T 

E 

X 

E*{MeV) 

0 

0 

0 

-112 

0 

2 

0 

0 

-94 

3.36 

4 

0 

0 

-52 

11.19 


(b) [f] = [3,l] (A,p) = (2,l) -contains the scissors mode (L = 1, S' = 0, 
T = 1). 

Note that the (S',T) combinations (0,1), (1,0) and (1,1) are allowed. 


L 

E 

X 

E*{MeV) 

1 

-58 

10.07 

2 

-46 

12.31 

3 

-28 

15.67 
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(c) [f] = [2,2] (A,/i) = (0,2) The {S,T) combinations (0,0), (0,2), (2,0) and 
( 1 , 1 ) are allowed. 


L 

E 

X 

E*{MeV) 

0 

-40 

13.43 

2 

-22 

16.78 


(d) [f] = [2,l,l] (A,/i) = (l,0) The {S,T) combinations (0,1), (1,0), (1,1), 
( 1 , 2 ) and ( 2 , 1 ) are allowed. 

L I E*{MeV) 

1 -10 19.02 


Note that this supermultiplet also has a state with the quantum numbers 
of the scissors mode L = 1S' = 0T = 1. 

Some further comments are in order. The scissors mode state in 
as a single band state originally discovered in electron scattering |^, was 


found when hner resolution ( 7 , 7 ') experiments were performed to consist of 
many states |T 8 |, . This was a beautiful example of intermediate structure. 

The supermultiplet scheme here affords a concrete example of the origin of 
the intermediate structure. Our scissors mode state at an energy of —58 x is 
degenerate with anL = lS'=lT = 1 state. If we introduce spin-dependent 
interactions the two states will admix and the degeneracy will be removed. 
We will get intermediate structure. 


8.2 B. Supermultiplet Scheme in 

We now give the energies and some properties of the states in 
(x = 0.1286): (a) [f] = [4,2] (A,/i) = (2,2) (includes ground state). 
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Allowed states: 


L 


T 

E 

X 

E*{MeV) 

0 

0 

1 

-96 

0 

2i 

0 

1 

-78 

2.32 

22 

0 

1 

-78 

2.32 

3 

0 

1 

-60 

4.63 

4 

0 

1 

-36 

7.72 


Note that the 2f and 2^ states are degenerate in this scheme. This arises 
from the fact that, in a rotational pictnre, the K valnes that are allowed are 
/i, /i — 2, .... Thns we have a degeneracy of J = 2^^ A' = 0 and J = 2^ 
K = 2. This degeneracy does not correspond to the experimental sitnation 
-the 2f and 2^ states are at 3.368 MeV and 5.960 MeV respectively -well 
separated. 

As a practical matter this degeneracy gives problems for the shell model 
code OX BASH |^. Shell model rontines often give wrong answers for 
transition rates when the states involved are degenerate. To overcome this 
difhcnlty, we have introdnced weak additional terms in the Hamiltonian to 
remove the degeneracy e.g. we nse a weak one-body spin-orbit interaction 
■ s with ^ = 0.1 MeV. 

(b) Two degenerate bands [f] = [4,l,l] (A,/i) = (3,0), [f] = [3,3] (A,p.) = (0,3) 

Allowed states: 


L 

A 

T 

E 

X 

E*{MeV) 

1 

1 

1 

-66 

3.853 

3 

1 

1 

-36 

7.716 
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We get the low-lying 1+ states 

(one 

from [4,1,1] from [3,3]). Note 

however that we 

have L = 1, S = 

= 1, hence the states cannot be excited by 

either the orbital operator or the spin operator. 


(c) Band which contains the scissors mode [f] = 

= [3,2,1] (A,/i)=(l,l) 

L 

S 

T 

E 

X 

E*{MeV) 

scissors mode 1 

0 

1 

-30 

8.49 

scissors mode 1 

0 

2 

// 

// 

1 

1 

1 

// 

// 

1 

0 

2 

// 

// 

1 

2 

1 

// 

// 


There are also several L = 2 states with ^ = —18 and E* = 10.03 MeV. 

X 

Note that the T = 1, S' = 0 scissors modes hnally make their appearance. 
There are two branches -one with isospin T = 1 and one with isospin T = 2. 
These two scissors mode states are degenerate in energy in the supermultiplet 
scheme ^^Be. 

(d) [f] = [2,2,2] (A, z/) = (0,0) 

L I E*{MeV) 

0 0 12.34 


The (S T) values are (0 1), (1 0), (1 2), (2 1), (0 4) and (4 0). 

8.3 B(M1) Transitions in the SU{3) Scheme 

For ®Be and the strength for orbital Ml transitions from J'^ = O’*" 
to scissors mode states can be obtained in the SU{3) scheme by observing 
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that in the process the ground state intrinsic state is transformed to the 
corresponding intrinsic state of state. For example in ®Be, the orbital 
isovector Ml operator, tI = {L'" — L'^) ,transforms the four nucleon 

intrinsic state (4,0) into the intrinsic state of I"*"! state that is (21). The 
operator (Lg — Tg) is the generator of the scissors mode. The orbital part 
of ground state and I'*'! state of ®Be in SU{3) scheme can be projected out 
from the corresponding maximum weight intrinsic states | [/](A,/i)e. A, p > 
where e = 2A + /i and A = p = A/i, by using the following projection]^, 

(2L + 1) 


\[f]{X,p)K,L,M>= 


dQD^^j,{Q)R{Q)\[f]{X,p)e,A,p> . (1) 


a{XpLK) 

Eq.(|l|) is a general equation for projecting out orbital part of the wave func¬ 
tion for a given L from a (A, p) state. In particular for ®Be we have, 

1 


and 


|[4](40)0,0,0>= 


|[31](21)1,1,0>= 


a(4000) 


dnD'^Q{n)R{n) \ [4] (4o)8, o, o > ( 2 ) 


<inD;,(Ji)fl(SJ)||31|(21)5,i 1 >. (3) 


a(2111) 

Knowing the [f] representation of the (A, p) states one can look up the cor¬ 
responding conjugate charge spin states to get dehnite JT states. 

The orbital magnetic transition strength B(M1) between these states is 
calculated by a method similar to that outlined in Appendix A of ref, pi 
and is found to be 


B{M1) = i|;;| < [31](21)110, ^ = 0,T = 1|(L^ - Lg")|[4](40)000, 5 = 0,T = 0 > 

= f An = 0-637 pI 

(4) 

In an analogous fashion one can calculate the scissors mode Ml transition 
strengths for the nucleus ^^Be. We hnd the following results for SU(3) limit 
transition strengths in ^^Be: 
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B(M1)(0+1 ^ 1+1) = = 0.0895^J 

B(M1)(0+1 ^ 1+2) = = 0.1492,,J. 

8.4 Realistic Spin-Orbit Interaction and Restoration 
of SU (3) Symmetry 

As pointed out before, an important role of spin dependent part of in¬ 
teraction is to remove the degeneracies present in the SU (3) limit by mixing 
up the same hnal angular momentum states arising due to a given intrinsic 
state as well as from different intrinsic states. In a realistic interaction, the 
relative strengths of spin dependent and spin independent part of interaction 
determine whether the wavefunctions are close to SU(3) scheme or a (j-j) 
coupling scheme is a better description of the system. 

To understand further the part played by spin independent part of the full 
realistic interaction in the restoration of SU (3) symmetry, we consider a small 
space calculation with the full spin-orbit part of the (x, y) interaction plus a 
variable Q.Q interaction. Figure.(7) is a plot of isovector orbital, spin and 
total strength for Ml transitions from J = = 0 —J = l^T = 1 states 

versus f, the parameter multiplying the full Q.Q interaction matrix elements 
for ®Be. For the spin part we use the operator 9.412X] <ytz i.e. we include the 
large isovector factor. In ®Be, with increasing t the orbital isovector strength 
is seen to approach the SU(3) limit value of 0.637p^. The contribution of 
isovector spin transition, on the other hand to total B(M1) decreases as t 
becomes large. This is because the SU{4:) limit is being approached and in 
this limit the spin contribution vanishes. 

In Fig. (7) with only spin-orbit part of realistic interaction in play, the cal¬ 
culated Ml transition strength for 0’''0 to I"*"! transitions has a large spin flip 
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contribution and is found tobe as large as 9.7 In a small space calculation 
with full realistic interaction(a;, y)(Table I), a total B(M1) value of 1.0547yU^ 
is obtained with an isovector orbital transition strength of 0.67/i^ and an 
isovector spin contribution of 0.38/i^. Of course to get the physical B{M1) 
we add the spin and orbital amplitudes and square. The spin B{M1) is a 
factor of 25 lower here than the t = 0 value in Fig 7. It still has some effect 
because of the factor 9.412. We may note that the orbital transition strength 
arising due to full realistic interaction is very close to the SU{3) limit indi¬ 
cating that the realistic interaction favors a restoration of SU (3) symmetry. 
The large space realistic interaction calculation inspite of the correlations 
induced by shell mixings results in a B(M1) value 1.2866/i^ and isovector 
orbital transition strength of 0.728/i^ indicating that the wavefunctions are 
still very close to SU (3) wave functions. 

In ^°Be the situation is more interesting due to splitting of scissors mode 
strength into two degenerate states in SU(3) limit. Figures (8) and (9) show 
the orbital and spin part respectively of Ml transition strength for transitions 
from ground state to J = l’''T = 1 states, J = 1~^T = 2 states and all J = 
states. For ground state to J = l’''T = 1 transitions the orbital B(M1) is seen 
to dip to a minimum for t = 0.3 before it starts increasing so as to approach 
its SU(3) limit. An opposite trend is observed in the corresponding spin 
strength that shows some increase, reaches a maximum and then tends to 
the SU{4) limit of zero. The characteristic behaviour at t = 0.3 is possibly 
a manifestation of a shape change at small deformation before the system 
stabilizes by acquiring a permanent deformation. The Ml transition sums 
for ground state to J = l’''T = 2 states, show a behaviour similar to that 
observed for ground state to J = l'''T = 1 transitions in ®Be. 
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9 Magnetic Dipole Transitions To Individual 
States 

We here present several tables of magnetic dipole transitions from the 
J = 0+ ground states of ^Be and to individual J = 1+ states. We use 
both the {x,y) interaction with x = 1, y = 1 and the Q ■ Q interaction. 

Concerning the latter, we learned in the previous section that there are 
many degenerate states in the Oku calculation when a Q ■ Q interaction is 
used. Unfortunately, most shell model routines, including the one used here, 
give erroneous results for transition rates when there are degeneracies. In all 
our small space (Oho;) calculations using Q ■ Q we have added a small spin- 
orbit interaction —■ s with ^ = 0.1 MeV. This works but it introduces 
an artihcial complexity in our tables. However, it is easy to see by eye 
what states would be degenerate if the spin-orbit interaction is removed. 
Alternatively, one can use the analytic expressions for the energies in the 
previous section. 

The columns in Tables VI through XIII are dehned as in Table I: 


(a) Isovector Orbital 

(b) Isovector Spin 

(c) Physical 


QItt 9lu 9 stt 9su 

0.5 -0.5 0 0 

0 0 0.5 -0.5 

1 0 5.586 -3.826 


9.1 Calculated Magnetic Dipole Transitions in ^Be 

In Tables VI and VII we present the details of the Qhoj calculated B{M1) 
values from the J = 0^, T = 0 ground state of ^Be to the J = 1^, T = 1 
excited states. 
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For the realistic {x, y) interaction, the isovector orbital strength is con¬ 
centrated in three states at 13.7, 16.6 and 18.0 MeV. The sum of the orbital 
strengths is |= 0.67 This is in fair agreement with the ex¬ 

perimental value B{M1) |= 0.81 (which is actually deduced from the 
downward 7 decay of the 17.64 MeV J = T = 1 state to the ground 
state). However, in the experiment all the strength is concentrated in one 
state whereas in our calculation we have considerable fragmentation. On 
the other hand, if we look at the physical transitions, there is much more 
concentration in one state at 13.73 MeV with B[M1) ]= 0.72 /i^. We will 
discuss this more soon. 

With the Q ■ Q interaction, all the orbital strength is concentrated in the 
(2-fold degenerate) state at 10.1 MeV with a summed strength B{M1) |= 
0.64 /i^, very similar to that of the {x,y) interaction. The energy is too 
low compared with experiment, but we must remember that we did not 
renormalize the strength y to allow for AiV = 2 admixtures. Note that the 
isovector spin transitions are zero with the Q ■ Q interaction. This is because 
we are at the SU (4) limit since Q ■ Q is a spin-independent interaction. 

Note that with the Q-Q interaction the summed orbital strength is ^ 
conhrming the expressions that were derived in the previous section. Coming 
back to the {x, y) interaction, we see that here also the isovector spin tran¬ 
sitions are very weak. But note that for the 13.73 MeV state whereas the 
orbital value of B{M1) ] is 0.2569 and the spin value is 0.0013 /i^, the 
physical value is 0.7155 jji\. The reason is that the spin and orbit amplitudes 
add coherently and that the spin amplitude is multiplied by the factor 9.412. 

For other states there is destructive interference between spin and orbit. 
For example, for the 16.64 MeV state, the value of B{M1) \ is 0.234 jj,% for 
the orbital case but only 0.065 for the physical case. 
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9.2 Calculated Magnetic Dipole Transitions in 

In Tables VIII and IX we present the details of the Oku calculated 
B{M1) values from the J = O’*" T = 1 ground state of to the J = 
I"*" T = 1 and to J = I’*' T = 2 excited states. We caution the reader that 
whereas for ^Be we presented the results in units of (Tables VI and VII), 
for ^^Be we use 10“2 as the unit. The reason for this is that the orbital 
transitions to individual states in ^^Be are considerably smaller than those 
in ^Be. 

Let us look at the Q ■ Q interaction (Table IX) first. There are several 
outstanding features which are explained in the previous section on L S' 
coupling and supermultiplet symmetry. 

The hrst two J = 1+ T = 1 states are degenerate at E* = 3.86 MeV. 
They carry no spin or orbital strength from the from the ground state. The 
[f] symmetries are [4 1 1] and [3 3]. They have additional quantum numbers 
L = 1S'=1T = 1. Since the isovector orbital operator {L^-Li^) cannot 
change both L and S from zero to one, the orbital B{M1) vanishes. A 
similar argument holds for the isovector spin operator. These lowest two 
states are therefore not scissors mode states. 

Then we have a four fold set of degenerate states with [3 2 1] symmetry 
at about 8.5 MeV excitation which does include the L = 1S' = 0T = 1 
scissors mode. We note that for ^^Be, the T = 1 scissors mode is degenerate 
with the T = 2 scissors mode also at 8.5 MeV excitation. This again is a 
prediction of the supermultiplet theory. 

The summed isovector orbital strength is ^ from J = O"*" T = 1 to 
the J = 1+ T = 1 states, and it is ^ to the J = 1+ T = 2 states. We 
have in effect a (2T + 1) rule: 

(2T + 1)t=2 5 

(2T + 1)t=i ^ 3 
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This coincides with the ratio of T = 2 to T = 1 strength. 

Recalling that the ^Be strength was ^ we see that the ratio of total 
strength ^ is 

We now come back to Table VIII which shows the same calculational 
results with the ‘realistic’ (x, y) interaction. There are several similarities but 
also some differences with the Q-Q results. Just as with the Q-Q interaction, 
the orbital transitions to the lowest two J = !■*■ T = 1 states at 6.14 and 7.68 
MeV are very weak 0.16 and 0.17 x (10“^/r^) respectively. However, the 
spin transitions, which with Q ■ Q were also zero, are now sufficiently strong 
so as to have a visible effect. For example, the physical B{M1) | to the 7.68 
MeV state is calculated to be 1.85 This is certainly measurable. 

As with the Q ■ Q interaction, the scissors mode states with the (x, y) 
interaction are at a much higher energy than the lowest two 1^ states 19 
MeV. Also, the J = 1+ T = 1 and T = 2 excitations are roughly in the 
same energy range -the Q ■ Q interaction has them degenerate. The ratio of 
T = 2 to T = 1 orbital strength is about the same for the {x, y) interaction 
as ioT Q ■ Q 1.44 vs. |. 

One major difference is that the energy scale is larger for the {x, y) inter¬ 
action than for Q ■ Q. The lowest and higher energies in Table VIII are 6.14 
MeV and 30.96 MeV whereas in Table IX they are 3.85 MeV and 12.35 
MeV. 

In Tables X and XI we present results of large space calculations for 
to be compared with the corresponding small space Tables VI and VII. 
Likewise in Tables XII and XIII we present the large space results for ^^Be 
to be compared with tables VIII and IX. We do not show all the states here, 
only the low energy sector. The excitation energies are in general larger in 
the large space calculations. The major changes occur when one has nearly 
degenerate levels sharing some strength. For example, the lowest two 



states in the large space calculations, which are at 7.38 MeV and 9.62 MeV 
have almost equal Ml strengths 0.64 and 0.79 In the small space, 
the lowest state has only 0.011 and the second one 1.8 

A sensible attitude is to assume that neither calculation is accurate enough 
to give the detailed distribution of strength between the two states -only the 
summed strength for the two states should be compared with experiment. 

With the Q ■ Q interaction in the large space, the degeneracy encountered 
in the small space calculation is removed. In part, this is due to the fact 
that we do not include the single-particle terms J2i=j Q{i) ■ Q{j)- These will 
induce a single-particle splitting between Is and Od in the N = 2 shell. But 
since degeneracies give us trouble in our shell model diagonalizations, we are 
happy to leave the calculation as is. 

With Q-Q the scissors mode strength in ^Be gets pushed up from the small 
space value of 10.1 MeV to the large space value of 15.5 MeV. For ^^Be the 
corresponding numbers are 8.5 MeV and 11.3 MeV for the J = 1^, T = 1 
states and essentially the same for J = I’*', T = 2 states. That is, the near 
degeneracy of the T = 1 and T = 2 scissors modes in ^^Be is maintained in 
the large space calculation. 

Lastly, we reiterate the fact that the shell model calculations here yield 
not only colective magnetic states but also show intermediate structure. That 
this structure is a natural occurence is shown by the supermultiplet model, 
where for a given [/]l=i there are several S and T values possible. It is of 
course very difficult to get the details of the intermediate structure to come 
out right, but it is good to be able to explain the origin of this structure. 
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10 Additional Comments and Closing Remarks 


We can gain fnrther insight into the natnre of ^Be and by evaluating 
the quadrupole moments of the J = 2+ states. A small space calculation gives 
the following values: 


^Be J = 2+, T = 0 
lOfie J = 2t, T = . 
^^Be J = 2t, T = : 


Q ■ Q — 0.1/ ■ s 
Q = —8.02 e fm 
Q = —2.52 e fm' 
Q = +2.06 e fm‘- 


{x, y) interaction 
Q = —7.86 e /m‘ 
Q = —7.68 e /m‘ 
Q = +6.91 e fm 


In the rotational model the quadrupole moment of the 2+ of a 77 = 0 band 
is —|Qo where Qq is the intrinsic quadrupole moment. Thus a negative Q 
corresponds to a prolate shape and a positive Q to an oblate shape. From 
the above, ^Be acts as a normal deformed nucleus of the prolate shape. 

It has been pointed out by Harvey that in the SU{3) scheme, whenever 
fi is less than or equal to A the nucleus becomes oblate [|l^. For the ground 
state band in ^Be A is bigger than y but for ^^Be A and y are equal. The 
situation with ^^Be is somewhat confusing. With the Q-Q interaction, which 
one might think would favor deformation, the quadrupole moment of the 2f 
state drops to -2.52 e fmf. Recall that for a perfect vibrator, the value of Q 
is zero, so it would appear that ^^Be is headed in that direction. However 
with the realistic interaction, which contains a large spin-orbit interaction 
that one might think would oppose deformation, the quadrupole moment of 
^^Be becomes more negative -almost the same as that of ^Be. Note also 
that the calculated values of Q for the 2i and 2f states are nearly equal but 
opposite to both interactions. 

We have learned many interesting things by considering scissors modes in 
light nuclei. First of all there is evidence for their existence. This evidence 
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comes strangely from a nucleus whose ground state is unstable -^Be. We learn 
of the existence from the inverse process i.e. 7 decay from the J = 1^, T = 1 
state at 17.64 MeV to the ground state The decay to the 2i state, 
presumably a member of a 77 = 0 rotational band, is also observed and this 
suggests that theoretical studies (and experimental ones as well whenever 
possible) should be made not only between between J = 1 + 77 = 1 and 
J = O’*" 77 = 0 states but also between J = 1+ 77 = 1 and J = 2+ 77 = 0 
states. This will make the picture of scissors modes more complete. In this 
work we considered but one example and showed that the ratio of J = I’*' 
decay rate to J = 2 + vs J = 0 ^ deviates from the simple rotational formula 
result of 0.5. Further studies along these lines are being planned. 

To make the picture even more complete, one can also study the decay 
of J = I"*" 77 = 1 to J = 2+ 77 = 2. We were almost forced into such a 
study by the fact that in the SU{3) scheme there is a two-fold degeneracy of 
the lowest J = 2 + states in [0 0 - Presumably, these two states are 
admixtures of 2+ 77 = 0 and 2+ 77 = 2. 

The shell model approach used here 113 enables us to study fragmen¬ 
tation or intermediate structure. We hnd for example that with an elec¬ 
tromagnetic probe there are, besides the strong scissors mode states, some 
almost ‘invisible’ states. These have separately substantial orbital contribu¬ 
tions and substantial spin contributions to the magnetic dipole excitations 
but the physical B{M1) is very small because of the destructive interference 
of the spin and orbital amplitudes. For example, as seen in Table X, in ®7?e 
we calculate that the low lying orbital strength is shared almost equally be¬ 
tween two states (at 18.0 MeV and 20.8 MeV) -the strengths being 0.22 

and 0.30 respectively. However, the physical 7?(Ml)’s are 0.62 
and 0.12 /i|f. Thus one can miss considerable orbital strength into these 
‘invisible’ sttaes if one uses only an electromagnetic probe. 
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Another thing we learn is that although spin excitations are strongly 
suppressed they cannot be ignored. In the SU{4:) limit, the spin matrix 
elements vanish and there is a clear tendency in our calculations in that 
direction. However, since the isovector spin operator is multiplied by a factor 
of 9.412, the spin and orbital contributions tend to be on the same footing. 

In the example of the above paragraph in the decay of the (calculated) 
18.0 MeV state, the orbital B[M1) is only 0.22 but the physical one which 
induces the spin contribution is 0.62 /i^. In the hrst two J = 1+ states 

are calculated to be excited mainly by the spin operator and the i?(Ml)’s 
should be substantial 0.5 /i^. Yet in the U{?>) — S'f/(4) limit, these lowest 
two states [f] = [4,l,l] and [3,3] should not be excited at all either by the spin 
or by the orbital operators. 

We have found the Wigner supermultiplet scheme combined with the 
SU (3) scheme |^, a very useful guide to the more complicated shell model 
calculations that we have performed. There is the added simplicity in the p 
shell that there is a one-to-one correpondence between a given [f] symmetry 
and the (A, p) symmetries. Many interesting properties about scissors modes 
can be literally read off the pages of the book by Hammermesh [113. For 


example, there is the fact that the T = 1 and T = 2 scissors mode states 
in ^^Be are degenerate in energy. This is an exact result with the Q ■ Q 
interaction in a p shell calculation. Results very close to this are obtained 
with a realistic interaction, but we frankly didn’t notice this until we made 
an SU{3) analysis. Also the non-obvious fact that the lowest two I"*" states 


m 
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Be are not scissors mode states is made clear by such an analysis. 

Also the fact that scissors mode states everywhere, including have 


intermediate structure [^, [1^ is made clear by the supermultiplet scheme. 
For a given L = 1 state there are often several S, T combinations which are 
degenerate. The removal of the degeneracy and the mixing of these states e.g. 
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by a spin-orbit interaction leads to fragmentation and intermediate structure. 

By extending the shell model calculations to ‘large space’ i.e. by including 
2huj excitations, we were able to calculate the cumulative energy weighted 
strength distribution for isovector orbital excitations. The results which are 
shown in several hgures are characterized by a low energy rise followed by a 
second plateau. The shapes of the distributions were similar for the two con¬ 
trasting interactions used here -the ‘realistic’ short range interaction and the 
schematic quadrupole-quadrupole interaction. The results were compared 
witht the simple Nilsson model ii which predicts that the energy-weighted 
sum at high energy (beyond the hrst plateau) should equal the low energy 
rise i.e. the ratio ■;—— should be 2:1. The actual calculated ratios witht 

LOW energy 

the {x,y) interaction are 1.75 for ^Be and 2.52 for ^^Be. The correspond¬ 
ing numbers for the Q ■ Q interaction were 1.37 and 2.33. We see that the 
Nilsson model is not bad as a hrst orientation but there are huctuations. 
The fact that the above ratios are larger for the (x, y) interaction than for 
the Q ■ Q interaction may support the idea of Hamamoto and Nazarewicz [§ 
that the symmetry energy will cause the ratio to increase. Our calculations 
however do not support their claim that the high energy part of the energy 
weighted orbital strength should always be much larger than the low energy 
part -certainly not for a ‘normal’ rotational nucleus like ®i?e. For ^^Be the 
ratio is however somewhat larger than the Nilsson model prediction. 

Whether this is due to the atypical properties of ^^Be mentioned in the 
text or is a harbinger of what will happen for most other nuclei remains to be 
seen. From an experimental point of view, it would be helpful to have more 
data on ^^Be. Not only have no J = !■*■ states been identihed but neither 
has the J = 4^ state been seen. The location of this state might help us 
decide whether ^^Be is rotational or vibrational. 
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At any rate, we should examine a larger range of nuclei and look into more 
detail about the symmetry energy in order to be able to make more dehnitive 
statements about the systematics of the cumulative energy weighted distri¬ 
butions. We note that the Zheng-Zamick sum rule is able to handle the 
divergent behaviour between ^Be and ^^Be. This sum rule involves the dif¬ 
ference between isoscalar and isovector summed B{E2) strength, whereas 
corresponding expressions by Heyde and de Coster p based on the LB.A. 
model ^ as well as empirical analyses involve only B{E2) to the 
lowest 2+ states. Even here more sharpening up is in order. 

Our initial reason for studying lighter nuclei is that they would give us 
insight into the behaviour of heavier nuclei, and we could carry out more 
complete calculations in the low A region. But then we found many results 
which made light nuclei studies fascinating for their own sake. One rather 
broad lesson we have learned in the light nucleus study is that there can 
be considerable change in going from one nucleus to the next, and perhaps 
in heavier nuclei too much effort has been made to make the nuclei £t into 
a smooth pattern. For example, we hud that whereas in ^Be the lowest 
state is dominantly a scissors mode state, in ^^Be the lowest states 
are not scissors mode states at all -they can only be reached by the spin 
operator. We should perhaps be looking for more variety of behaviour in 
heavier nuclei. Lastly the supermultiplet scheme which we found extremely 
useful was originally thought to be of interest only in light nuclei where the 
spin-orbit interaction is small relative to the residual interaction. However, it 
is now being realized that even in heavy nuclei -especially for superdeformed 
states this scheme may once again be very relevant. This would make our 
light nuclear studies all the more important. 
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Table I. Summed Magnetic Dipole Strengths (in 


The (x,y) interaction (x = 1, y = 1) 

Space 

Isovector Orhital^^^ 

Isovector Spin^^^ 

Isoscalar Orhital^^'> 

Physical^^^ 

^Be J = of T = 0 ^ J = 1+ T = l 


0.6701 

0.00427 


1.0547 

Large^^'^ 

0.7283 

0.00622 


1.2866 

Low Large^^^ 

0.5890 

0.00322 


0.8999 

i°Be J = 0+T = 1^J = 1+T=1 

Small 

0.1112 

0.0234 

0.0245 

2.0930 

Large 

0.1963 

0.0208 

0.0270 

1.9517 

Low Large 

0.1002 

0.0187 

0.0200 

1.6070 

i°Be J = 0fT = l^J = l+T = 2 

Small 

0.1508 

0.00105 


0.0597 

Large 

0.1830 

0.00222 


0.2276 

Low Large 

0.1339 

0.000928 


0.0754 

The Q ■ Q interaction 

^Be J = 0^ T = 0 ^ J = 1+ T = 1 

Small 

0.6364 

0.0000 


0.6364 

Large 

0.7408 

0.0005 


0.7858 

Low Large 

0.6593 

0.0002 


0.6775 

i°Be j = 05 ^r = i^j = i+r=i 

Small 

0.0895 

0.0001 


0.0986 

Large 

0.1788 

0.0004 


0.2044 

Low Large 

0.0922 

0.0000 


0.0881 

i°Be J = 05 ^r = l^J = l+r = 2 

Small 

0.1492 

0.0000 


0.1486 

Large 

0.1744 

0.0002 


0.1950 

Low Large 

0.1513 

0.0000 


0.1617 
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Table I Captions 

(a) Small Space (Os)^(Op)® 

(b) Large Space (Os)^(Op)® + all 2hoj excitations 

(c) Low energy part of Large Space (up to the first plateau) 

(d) gin = 0.5 giu = -0.5 gsn = 0 gsu = 0 
(®) 9itt 0 gii/ 0 ggn 0.5 ggi, 0.5 

(f) gi-K = 0.5 gii, = 0.5 gsn = 0 ggy = 0. The value for isoscalar spin is the same 
as for isoscalar orbital in the case of J = 0’'', T = l—T=l. 
For AT = 1 the isoscalar case gives zero. 

(g) gin = l giu = ^ gsn = 5.586 gg^ = -3.826 
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Table II. Summed Energy Weighted Magnetic Dipole Strengths (in nj^MeV) 


The {x,y) interaction (x = 1, y = 1) 

Space 

Isovector Orbital 

Isovector Spin 

Isoscalar Orbital 

Physical 

^Be J = 0t T = 0 ^ J = 1+ T = l 

Small 

10.689 

0.08359 


16.878 

Large 

20.744 

0.2789 


44.092 

Low Large 

11.854 

0.0778 


18.349 

i°Be J = 0+T=l^J = l+T=l 

Small 

2.0750 

0.1954 

0.2328 

18.852 

Large 

7.9511 

0.3119 

0.6929 

35.909 

Low Large 

2.4680 

0.1814 

0.2369 

18.429 

i°Be J = 05 ^r = l^J = l+r = 2 

Small 

2.984 

0.0229 


1.4783 

Large 

6.720 

0.1167 


12.83 

Low Large 

3.3590 

0.0267 


2.2843 

The Q ■ Q interaction 

®Be J = 0^ T = 0 ^ J = 1+ T = 1 

Small 

6.411 

0.0000 


6.411 

Large 

14.040 

0.0229 


15.861 

Low Large 

10.282 

0.0037 


10.619 

i°Be J = 0fT=l^J = l+T=l 

Small 

0.7597 

0.0004 


0.7942 

Large 

3.8108 

0.0138 


4.8627 

Low Large 

1.041 

0.0029 


0.9904 

i°Be J = 05 ^r=l^J = l+T = 2 

Small 

1.2661 

0.0000 


1.2619 

Large 

2.6019 

0.0068 


3.1687 

Low Large 

1.712 

0.0005 


1.8335 
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Table III. The B{E2) values in e^fm^ in ^Be with (a) The {x,y) interaction 
with X = 1, y = 1 and (b) the Q • Q interaction. 


Interaction 

mam 


o 

II 

T = 0^J = 2+T = 0 

2f Small 

17.59 

18.13 

2f Large 

25.97 

49.16 

Sum Small 

17.82 

18.13 

Sum Large 

35.21 

59.38 

Sum Low Large 

33.28 

49.22 

j = ot 

r = i^j = 2+r = i 

Sum Small 

2.31 

2.59 

Sum Large 

13.47 

22.44 

Sum Low Large 

2.10 

13.51 
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Table IV. The B{E2) values in e^fm^ in (large space) with (a) The (x, y) 
interaction with x = 1, y = 1 and (b) the Q ■ Q interaction. 


Interaction 

x = l,y = l 

^muo'^r'^ — xQ ■ Q 

J = 0^ 

T=1^J = 2+T = 1 

2f Small 

22.98 

17.31 

2^ Small 

0.30 

7.67 

Sum 

23.28 

24.98 

2f Large 

19.68 

5.17 

2^ Large 

3.22 

41.23 

Sum 

22.90 

46.40 

Sum Small 

24.21 

25.50 

Sum Large 

39.74 

74.70 

Sum Low Large 

24.55 

72.14 

o 

II 

r = i^j = 2+r = 2 

Sum Small 

0.7858 

0.8016 

Sum Large 

5.851 

12.19 

Sum Low Large 

0.5669 

10.57 


40 








Table V. Ratio of orbital Ml strength and orbital energy weighted Ml strength 
for the (x,y) interaction with x = \, y = 1. 


Space 

Ratio of orbital 

Ratio of energy weighted 


strength 

orbital strength 

Small 

0.3910 

0.4730 

Large 

0.5208 

0.7072 

Low Large 

0.3974 

0.4920 


Table VI. Calculated Magnetic Dipole Strength in ^Be with the (x, y) 
interaction (x = 1, y = 1) in small space (in units of fi%). 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 


= 0f T = l^ J 

^ = 1+ T=1 


13.73 

0.2569 

0.0013 

0.7155 

16.64 

0.2344 

0.0006 

0.0651 

18.05 

0.1748 

0.0002 

0.0778 

20.73 

0.0013 

0.0012 

0.1307 

26.80 

0.0007 

0.0005 

0.0319 

28.46 

0.0007 

0.0004 

0.0272 

29.66 

0.0014 

0.0000 

0.0001 

34.76 

0.0000 

0.0001 

0.0065 

SUM 

0.6701 

0.0043 

1.0547 
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Table VII. Calculated Magnetic Dipole Strength in ^Be with the Hamiltonian 
H = -^ + ^muP'r'^ — xQ ' Q with x = 0.5762 in small space 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 


II 

o 

II 

h-‘ 

i 

^ = 1+ T=1 


10.06 

0.4599 

0.0000 

0.4851 

10.11 

0.1765 

0.0000 

0.1524 

12.33 

0.0000 

0.0000 

0.0001 

13.43 

0.0000 

0.0000 

0.0000 

16.84 

0.0000 

0.0000 

0.0000 

18.97 

0.0000 

0.0000 

0.0000 

19.05 

0.0000 

0.0000 

0.0000 

19.11 

0.0000 

0.0000 

0.0000 

SUM 

0.6364“ 

0.0000 

0.6376 


(a) The summed isovector orbital strength is ^ /x^. 
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Table VIII. Calculated Magnetic Dipole Strength in with the (x, y) 
interaction (x = 1, y = 1) in small space (in units of 10“^ P^n)- 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 

A. J = 0| T = 1 ^ 

J = l+ T = 1 


6.14 

0.1641 

0.0578 

1.0740 

7.68 

0.1699 

2.1690 

184.90 

14.31 

0.1000 

0.0048 

1.2020 

16.48 

0.0000 

0.0105 

1.1050 

18.14 

7.2840 

0.0020 

10.530 

19.19 

1.7540 

0.0152 

4.8420 

22.62 

1.1930 

0.0064 

0.4793 

23.93 

0.3904 

0.0513 

2.6420 

25.72 

0.0606 

0.0231 

2.5720 

SUM 

11.120 

2.3409 

209.30 

B. J = 0^ T = 1 ^ 

J = 1+ T = 2 


17.91 

6.4260 

0.0331 

0.6757 

20.67 

7.1290 

0.0494 

0.3339 

23.50 

1.5280 

0.0039 

3.3170 

30.96 

0.0001 

0.0188 

1.6430 

SUM 

15.083 

0.1052 

5.9696 
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Table IX. Calculated Magnetic Dipole Strength in with the Hamiltonian 
~ 2 m + \rnuj‘^r‘^ — xQ ' Q with x = 0.3615 in small space (10“^ fi%). 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 

A. J = of T = 

1 ^ 

J = l+ T = 1 


3.848 

0.0006 


0.0051 

0.4173 

3.870 

0.0006 


0.0046 

0.5161 

8.444 

3.0350 


0.0001 

2.7460 

8.486 

0.0256 


0.0000 

0.0206 

8.491 

3.3300 


0.0000 

3.3190 

8.533 

2.5590 


0.0001 

2.8440 

9.988 

0.0002 


0.0000 

0.0001 

10.035 

0.0003 


0.0000 

0.0003 

10.082 

0.0001 


0.0000 

0.0002 

SUM 

8.9555“ 

0.0099 

9.8636 

II 

hn 

o 

II 

1 ^ 

J = 1+ T = 2 


8.441 

7.6680 


0.0003 

6.8270 

8.536 

7.2500 


0.0002 

8.0320 

10.034 

0.0000 


0.0000 

0.0000 

12.349 

0.0000 


0.0000 

0.0000 

SUM 

14.918“ 

0.0005 

14.859 


(a) The summed isovector orbital strength to T = 1 states is ^ and to T = 2 
states it is ^ 
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Table X. Calculated Magnetic Dipole Strength in ^Be with the {x, y) interaction 
(x = 1, y = 1) in large space^“^ (in nnits of n%). 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 


T 

o 

II 

Sh 

o 

II 

= 1+ T = 1 


17.96 

0.2208 

0.0011 

0.6160 

20.80 

0.2970 

0.0004 

0.1213 

22.89 

0.0670 

0.0001 

0.0213 

27.04 

0.0000 

0.0010 

0.0935 

33.51 

0.0007 

0.0002 

0.0081 

34.57 

0.0000 

0.0003 

0.0250 

35.73 

0.0007 

0.0000 

0.0000 

40.22 

0.0001 

0.0000 

0.0036 

46.07 

0.0003 

0.0000 

0.0110 

46.41 

0.0000 

0.0000 

0.0000 

48.46 

0.0060 

0.0001 

0.0001 

49.61 

0.0019 

0.0000 

0.0023 

49.92 

0.0012 

0.0000 

0.0004 

50.32 

0.0004 

0.0000 

0.0036 

SUM 

0.5986 

0.0033 

0.9063 


(a) Large Space; all Ohu> configurations plus 2?uv excitations. 
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Table XI. Calculated Magnetic Dipole Strength in with the Hamiltonian 
H = -^ + — xQ ' Q with x = 0.5762 in large space 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 


j = o+r = o^j 

= 1+ T = 1 


15.49 

0.6012 

0.0000 

0.5513 

15.70 

0.0552 

0.0001 

0.1189 

18.42 

0.0003 

0.0000 

0.0005 

22.69 

0.0000 

0.0000 

0.0020 

25.43 

0.0000 

0.0000 

0.0004 

27.82 

0.0000 

0.0000 

0.0014 

28.14 

0.0000 

0.0000 

0.0002 

28.24 

0.0000 

0.0000 

0.0003 

35.24 

0.0009 

0.0000 

0.0013 

37.88 

0.0003 

0.0000 

0.0000 

39.30 

0.0001 

0.0000 

0.0001 

39.51 

0.0003 

0.0000 

0.0003 

40.17 

0.0010 

0.0000 

0.0010 

SUM 

0.6593 

0.0002 

0.6775 
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Table XII. Calculated Magnetic Dipole Strength in with the (x, y) 
interaction (x = 1, y = 1) in large space*^“^ (in units of 10 “^ ^at)- 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 

A. J = 0| T = 1 ^ 

J = l+ T = 1 


7.38 

0.006 

0.909 

63.78 

9.62 

0.223 

0.847 

79.47 

18.29 

0.015 

0.000 

0.026 

20.65 

0.047 

0.010 

0.768 

22.73 

6.087 

0.004 

10.40 

23.99 

1.620 

0.013 

4.401 

27.15 

1.411 

0.004 

0.788 

28.90 

0.086 

0.047 

3.528 

30.86 

0.017 

0.031 

2.831 

37.20 

0.450 

0.001 

0.849 

39.96 

0.002 

0.003 

0.150 

40.39 

0.058 

0.002 

0.028 

SUM 

10.02 

1.870 

160.7 

B. J = 0^ T = 1 ^ 

J = 1+ T = 2 


23.27 

4.733 

0.015 

1.058 

25.99 

7.560 

0.046 

0.524 

29.11 

1.004 

0.013 

4.311 

38.33 

0.000 

0.017 

1.524 

49.42 

0.017 

0.000 

0.098 

49.67 

0.003 

0.000 

0.028 

50.84 

0.075 

0.000 

0.000 

SUM 

13.39 

0.093 

7.540 


(a) Large Space; all Oho; configurations plus 2fujj excitations. 
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Table XIII. Calculated Magnetic Dipole Strength in with the Hamiltonian 
H = ^ + ^muP'r'^ — xQ ' Q with x = 0.3615 in large space (10“^ 


^,(1+) MeV 

Isovector Orbital 

Isovector Spin 

Physical 


A. J 

= 0^ T = 1 ^ 

J = l+ T = 1 


3.60 


0.002 

0.000 

0.019 

6.45 


0.005 

0.001 

0.058 

11.19 


2.548 

0.000 

2.393 

11.31 


1.685 

0.000 

2.229 

11.32 


4.563 

0.001 

3.382 

11.35 


0.395 

0.001 

0.702 

12.92 


0.019 

0.000 

0.019 

12.98 


0.004 

0.000 

0.003 

13.02 


0.002 

0.000 

0.003 

22.69 


0.011 

0.001 

0.154 

23.12 


0.019 

0.000 

0.001 

25.42 


1.579 

0.000 

1.220 

25.94 


0.043 

0.003 

0.497 

26.51 


0.090 

0.000 

0.012 

26.65 


0.004 

0.000 

0.005 

SUM 

10.97 

0.007 

10.70 


B. J 

= 0| T = 1 ^ 

J = 1+ T = 2 


11.31 


13.82 

0.001 

15.79 

11.34 


1.301 

0.003 

0.348 

12.98 


0.000 

0.000 

0.002 

16.76 


0.000 

0.000 

0.007 

26.58 


0.002 

0.000 

0.001 

26.67 


0.002 

0.000 

0.010 

28.25 


0.001 

0.000 

0.003 

SUM 

15.13 

0.004 

16.16 

























Figure Captions 

Figure (1): The cumulative sum of the energy-weighted isovector orbital 
B{M1) strength for the Oi", 0 —!■*■, 1 transitions in ^Be with the realistic 
interaction (x = 1, y = 1). 

Figure (2): Same as Figure 1 but with the Q ■ Q interaction. 

Figure (3): The cumulative sum of the energy-weighted isovector orbital 
B{M1) strength for the Oi", 1 —> I"*", 2 transitions in ^^Be with the realistic 
interaction (x = 1, y = 1). 

Figure (4): Same as Figure 3 but with the Q ■ Q interaction. 

Figure (5): The cumulative sum of the energy-weighted isovector orbital 
B{M1) strength for the Oi", 1 —I"*", 1 transitions in ^^Be with the realistic 
interaction (x = 1, y = 1). 

Figure (6): Same as Figure 5 but with the Q ■ Q interaction. 

Figure (7): = 0|T = 0 —> J = 1+T = 1 states) versus t 

for ^Be. The solid line, dashed line and dot-dash line are the total, spin and 
orbital parts of . 

Figure (8): The orbital part of versus t for ^^Be. The solid 

line, dashed line and dot-dash line are the total, B{M1){J = OfT = 1 — 
J = 1~^T = 1 states) and Y. B{M1){J = OfT = 1 —J = 1+T = 2 states) 
respectively. 

Figure (9): Same as in Fig.(8) for spin part of YB{M1). 
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